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ABSTRACT 



A debris disk is a constituent of any planetary system surrounding a main sequence star. We study whether close stellar encounters 
can disrupt and strip a debris disk of its planetesimals in the expanding open cluster of its birth with a decreasing star number density 
over 100 Myrs. Such stripping would affect the dust production and hence detectability of the disk. We tabulated the fractions of 
planetesimals stripped off during stellar flybys of miss distances between 100 and 1000 AU and for several mass ratios of the central 
to passing stars. We then estimated the numbers of close stellar encounters over the lifetime of several expanding open clusters 
characterized by their initial star densities. We found that a standard disk, with inner and outer radii of 40 and 100 AU, suffers no 
loss of planetesimals over 100 Myrs around a star bom in a common embedded cluster with star density < 1000 pc^^. In contrast, we 
found that such a disk is severely depleted of its planetesimals over this timescale around a star born in an Orion-type cluster where 
the star density is > 20 000 pc^^. In this environment, a disk loses > 97% of its planetesimals around an M-dwarf, > 63% around a 
solar-type star, and > 42% around an A-dwarf, over 100 Myrs. We roughly estimate that two-thirds of the stars may be born in such 
high star density clusters. This might explain in part why fewer debris disks are observed around lower mass stars. 

Key words. Stars : circumstellar matter ; surveys ; stars: low-mass ; planetary systems : formation 



1. Introduction 

A debris disk surrounding a main sequence star is the collection 
of planetesimals (comets or asteroids) that are the leftovers from 
an early phase of planet formation. A debris disk is a constituent 
, of any planetary system in the core-accretion theory of planet 
. formation. It is the analogue of the Kuiper Belt beyond Neptune, 
' or of the asteroid belt between Mars and Jupiter, in our Solar 
System. The formation and evolution of debris disks and plane- 
tary systems are inter-related in theory, and can provide comple- 
mentary insight into each other. 

The dust content of a debris disk is the key to its observabil- 
ity. Small dust grains scatter the light of the central star and make 
I the disk visible at optical wavelengths. Large grains are heated 
' by the central star and efficiently re-radiate in the far-IR and 
(sub)millimeter domains where they are observed. It is thought 
that dust must be continually or episodicly replenished by mu- 
tual collisions between planetesimals because dust grains are re- 
moved on short timescales (a few millions years or less). Disks 
have to be several-fold dustier than the Kuiper belt (xlOO) to be 
detectable with our instrumental sensitivities, even for the near- 
est stars. Nonetheless, bright disks are detected around 32 + 5% 
of A-type dwarfs (Su et al 2006, Wyatt et al 2003), 16 ± 2.8% of 
solar-type dwarfs (Bryden et al 2006, Trilling et al 2008, Najita 
& Williams 2005) and <5% of M-type dwarfs (Lestrade et al 
2006 and 2009, Gautier et al 2007) according to searches for 
cold dust at A - IQfim by Spitzer and at (sub)millimeter wave- 
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lengths by radiotelescopes. Debris disks were recently reviewed 
by Wyatt (2008). 

The mechanisms responsible for enhanced collisional activ- 
ity and for grinding planetesimals to dust in debris disks remain 
unclear. The statistics of disk detections just recalled cannot be 
accurately interpreted unless these mechanisms are identified. It 
has been proposed that unseen planets dynamically stir the disk 
and generate copious amount of dust that evolve under gravity, 
radiation pressure, P-R, and stellar-wind drags (Moro-Martin & 
Malhotra, 2002; Wyatt 2003; Mustill & Wyatt 2009; Kennedy & 
Wyatt 2010). Steady-state equilibrium between collisional cas- 
cades of planetesimals and dust removal processes can explain 
how disks fade with age (Dominik & Decin 2003). Orbit res- 
onance crossing of giant planets can trigger abrupt showers of 
planetesimals in the system, severely depleting the disk and lim- 
iting its detectability (Gomes et al 2005, Morbidelli et al 2005, 
Tsiganis et al 2005). However, the distribution of disk fractional 
luminosities with star ages is inconsistent with the idea that this 
mechanism is common among debris disks (Booth et al. 2010). 
Finally, distant icy-planets can successively form in waves out- 
wards in the disk and generate dust rings (Kenyon & Bromley, 
2002). 

In this paper, we study whether stellar flybys during the first 
100 Myrs of the lifetime of a star, while it is still in the expand- 
ing open cluster of its birth, can significantly deplete a disk of 
its planetesimals, affecting its dust production and therefore de- 
tectability. Stellar flybys have been invoked to explain the high 
eccentricity orbits of some Kuiper Belt Objects such as Sedna 
(Kenyon & Bromley, 2004), the dynamics of planetary systems 
(Malmberg et al. 2007, Spurzem et al. 2009), and the structures 
of debris disks (Larwood, 1997; Mouillet et al. 1997; Kalas et 
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Fig. 1. Stellar flyby. The debris disk is non-self-gravitating, and has inner and outer radii of 40 and 100 AU at the start of the 
simulation. The central star is at the origin of the coordinate system. The trajectory of the passing star is parabolic, coplanar and 
prograde with respect to the disk. The mass ratio of the two stars is unity. At the closest approach, the miss distance is 200 AU and 
the resulting maximum velocity is 4.2 km/s. The fraction of planetesimals stripped off the disk and lost by the system is 13% during 
this close stellar encounter. 



al 2000; Kobayashi & Ida, 2001). Impacts of stellar flybys dur- 
ing the first few million years have also been studied in the con- 
text of the evolution of protoplanetary disks and planet formation 
(Bonnell et al 2001, Olczack, Pfalzner & Eckait 2009). 

For our study, we did not resort to a full N-body integra- 
tion of an expanding cluster with stars surrounded by disks of 
planetesimals. We instead divided the problem into parts by first 
estimating the fraction of planetesimals stripped off a disk by the 
passage of a star moving on a parabolic trajectory in Sect. 2, and 
second by estimating the probability of close stellar encounters 
in an expanding cluster in Sect. 3. Finally, we discuss our results 
in Sect. 4. 

2. Stripping a debris disk by stellar flybys 

In our model, the central and the passing stars are considered 
to be point masses, named iric and nip, respectively. The disk is 



made of 10^ massless planetesimals randomly distributed, radi- 
ally and azimutally, with the planetesimal surface number den- 
sity S(fl) = So(fl/ao) ''^^ Hence, the disk is not self-gravitating, 
i.e. the orbit of each planetesimal is computed under the sole 
gravitational attractions of m^. and nip in the restricted three-body 
problem. 

In our simulations, we considered two sizes for the disk : 
a standard configuration with inner and outer radii of 40 and 
100 AU, and a more compact configuration with respective radii 
10 and 40 AU, at the start of the computation. These char- 
acteristic radii are based on limited observations at present. 
Protoplanetary transition disks have inner radii ranging from 15 
to 73 AU and outer radii from 30 to 135 AU, as measured by the 
SMA (Andrews et al. 2011). More mature debris disks around 
main sequence stars have radii determined by means of either 
imaging or the fitting of their SED, and comprised between ~ 10 
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AU and ~ 300 AU (e.g. Dent et al. 2000, Greaves et al. 2005, 
Liseau et al. 2008, Corder et al. 2009, Marshall et al. 2011). In 
addition, disk size may depend on the central star mass; for ex- 
ample, the disk surrounding the M-star AU Mic is about three 
times smaller than the disk around the A-star fi Pic of the same 
age (12 Myr) (Augereau & Beust 2006, Augereau et al. 2001). 
However, a correlation between disk size and star mass has not 
yet been established because of a lack of sufficient data. Current 
surveys with Herschel in the far-IR (Matthews et al., 2010, Eiroa 
et al., 2010) will provide new and important constraints. We 
adopted the two disk configurations, standard and compact, in 
our simulations to help us prepare for any revisions of disk sizes 
in the future. 

Our first suite of simulations was conducted in a three- 
dimensional (3D) coordinate system where the disk is confined 
to the x-y plane but the parabolic orbit of the passing star can 
be in an inclined plane. The starting position of the passing star 
is chosen to be at a large distance (> 1000 AU) to ensure that 
the disk is not perturbed by an impulsive force at the start of 
the computation. The initial velocity for the passing star is set 
so that, at the closest approach, it reaches a given miss distance 
while on a parabolic trajectory. The miss distances were set to 
be between 100 AU and 1000 AU, and the resulting maximum 
velocities of the passing star at closest approaches were between 
3 and 12 km/s. This maximum velocity is related to the miss 
distance because of the parabolic nature of the orbit. 

The equation of motion of a massless particle under the grav- 
itational attractions of the central and passing stars were written 
in a coordinate system where the central star is at the origin. 
These equations were integrated numerically using the fourth or- 
der Runge-Kutta algorithm described in Press et al (1992). The 
numerical accuracy of our implementation of this algorithm was 
tested by computing the total energy of a test star, i.e. one not 
gravitationally interacting with the disk, set on a high eccentric- 
ity orbit with a small semi-major axis to produce the highest 
acceleration in the system. We set the time step for our integra- 
tion so that the total energy of this test star was conserved to an 
accuracy lower than 0.1 %. We also extensively compared our 
Runge-Kutta results to the ones that we obtained first with the 
computationally expensive Euler method. 

The dynamical evolution of our system is simulated over a 
few thousand years, long enough for the passing star to reach the 
closest approach and move away to a large distance. Our simula- 
tion is exemplified in Figure[T]for prograde and coplanar motions 
of the disk and passing star This example shows that the disk is 
first slowly perturbed and then significantly disrupted at closest 
approach where the tidal forces reach a maximum. Finally the 
disk settles in to its new configuration with bound planetesimals 
on elliptical orbits, while other planetesimals are either captured 
by the passing star or ejected on unbound orbits. We checked 
that the fractions of stripped planetesimals (captured -i- ejected) 
are accurately estimated for a disk confined to the x-y plane and 
accurately represented by 10"* particles in comparing to simula- 
tions with larger numbers of particles. 

Our simulations were carried out by varying the miss dis- 
tance from 100 AU to 1000 AU, the sense of circulation be- 
tween the disk and the passing star (prograde, retrograde), the 
inclination / of the orbital plane of the passing star (0°, 30°, 45°, 
90°), and the mass radtio niplnic (20, 10, 4, 2, 1, 0.5, 0.25, 0.1). 
We found that all stellar flybys with retrograde motions are in- 
efficient leading to very low fractions of stripped planetesimals. 
This is because they do not move in concert with the passing 
star and are not gravitationally perturbed for a sufficiently long 
period of time, as also noted by other authors in similar circum- 
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Fig. 2. Fractions of planetesimals stripped off a disk during stel- 
lar flybys at various miss distances and for various color-coded 
mass ratios mp/nic of the passing to central stars. The motions 
of the disk and of the passing star are prograde and their or- 
bital planes are coplanar The inner and outer radii of the disk 
are 40 and 100 AU, respectively, at the start of the simulation. 
Thick lines connect the points computed for a disk with plan- 
etesimals initially set on circular orbits. Dotted lines connect the 
points computed for a dynamically excited disk with orbits of 
planetesimals initially set with random eccentricities and incli- 
nations (see details at end of Sect. 2). 
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Fig. 3. Same as Figure |2]but for a compact disk with inner and 
outer radii of 10 and 40 AU. 



stances (Toomre & Toomre, 1972, Beust & Dutrey, 2006, Reche 
et al., 2009). All flybys with prograde motion and inclination 
/ > 45° are also inefficient. Only flybys with prograde motion 
and inclination / < 45° lead to significant stripping whose mag- 
nitude is quantified by Figures |2] and [3] These figures show the 
fraction of stripped planetesimals as a function of the miss dis- 
tance, for the mass ratios mp/iric mentioned above, for the stan- 
dard and compact disks, and for inclination / < 45° of the orbital 
plane of the passing star. Clearly, by comparing Figures |2] and 
|3] a compact disk with inner and outer radii of 10 and 40 AU is 
much less severely depleted of its planetesimals than the more 
extended one. 
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Table 1. Stellar flybys of miss distance 200 AU for two initial 
disk excitations *. 



nip/m^ 


Fractions of stripped planetesimals 




no excitation 


excitation 


4 


47% 


63% 


2 


27% 


42% 


1 


12% 


21 % 


0.5 


5% 


12% 


0.25 


1.2% 


6.2% 



' standard size disk. 



To check our simulations, we duplicated the results of 
Larwood and Kalas (2001) for the debris disk around/? Pic per- 
turbed by a passing star We found excellent agreement with their 
Table 1 for the numbers of particles unbound and captured by the 
passing star, for all ratios irip/mc and inclinations. We also com- 
pared our results with the study of stripping effects in the Kuiper 
Belt when the Solar System was in its birth aggregate performed 
by Adams and Laughlin (2001). They estimate that 40% of the 
Kuiper Belt objets with semi-major axis between 30 and 70 AU 
are removed from the Solar System in an encounter with a 1 so- 
lar mass star at the miss distance ~ 200 AU ~ 350/ V^ AU hy 
using their fiducial cross-section cr ~ (350 AU)^. This fraction 
is higher than our prediction of 13% in Figure |2] We ascribe this 
difference to their model which also includes the scattering of 
high eccentricity planetesimals by Neptune that is not consid- 
ered in our model. 

Finally, we studied how the initial conditions in the disk af- 
fect the outcome of a flyby since the disk can be initially excited 
by the dynamical perturbations of a previous flyby. We run a new 
suite of simulations in which planetesimals were set on eccen- 
tric and inclined, instead of circular, orbits. Their eccentricities, 
inclinations, and semi-major axes a were randomly distributed 
between 0.0 and 0.95, ±30°, and 40 and 100 AU, respectively, at 
the start of the simulation. The planetesimal surface number den- 
sity 2(a) was the same power-law as given above and lO*" par- 
ticles were used to represent this three-dimensional (3D) disk. 
These new simulations have yielded fractions of stripped plan- 
etesimals that are somewhat larger than for a disk with no initial 
excitation, but not exceedingly larger. These new fractions are 
the points connected by dotted lines in Figures |2] and [3] To ap- 
preciate the differences, we also present a subset of these points 
in Table[T]for a standard disk, miss distance of 200 AU and mass 
ratios nip/mc - 4, 2, 1, 0.5, and 0.25, which are representative 
examples. We estimated that the final depletions after 100 Myr 
between the two models can change by at most a factor of two 
depending on the exact history of encounters. Such a change is 
not important for the assessment of the impact of depletion on 
the dust production of a disk; only a depletion factor greater than 
ten is meaningful for the disk detectability. For the remaining 
analysis, we adopted the fractions given in Figures |2]and[3]for a 
disk with no initial excitation. 



3. Close stellar encounters in an expanding open 
cluster 

3. 1 . Number of stellar encounters over the cluster lifetime 

We do not resort to an N-body calculation to estimate the fre- 
quencies of encounters between stars in an open cluster. We 



instead use the kinetic theory complemented by gravitational 
focusing. In this framework, the encounter time tg„c can be 
estimated analytically by rederiving Eqs. 8-122 of Binney & 
Tremaine (1987, p. 539), which is valid for two stars of iden- 
tical masses, for the more general case of a central star and a 
passing star of masses m^ and nip moving in a field with uniform 
star number density n. The corresponding encounter rate at miss 
distance d < dene is 



— 4 y/n n cr de„ 



2 yJnGinic + nip) n d^. 



(1), 



where cr is the dispersion of the Maxwellian distribution in the 
stellar velocities of the cluster, and G is the gravitational con- 
stant. The first term is the collision rate inferred from the sim- 
ple kinetic theory of a star moving at the mean velocity 4cr/ a/tt 
(mean velocity of a Maxwellian with dispersion cr). This expres- 
sion is based on the volume swept by the cross-section Trii^,,^ at 
this mean velocity during te„c. The second term is gravitational 
focusing to account for the true orbits of the two stars. In conve- 
nient units, the total encounter rate is 



-iyr-^) = 1.9x10"*^' 



/ n \ / O- \ / dene V 

1 1000 pc-3 1 1 1 km/s / 1 100 AU / 



8.8xlO-''f^^^i^^]f " ]f^^)fi^^) (2). 
\ IMo AlOOOpc-VVlOOAUA cr I 

For our model, we assume that the star number density n is a 
linear function of time t. In our notation, this is 



■' nofp - (no - ni)fp(t/tci) 



(3), 



where fp is the stellar fraction for each spectral type range listed 
in Table |2] With this linear expression of time, the total star 
number density in the cluster decreases from no at f = to 
ni =0.1 pc~^ at ? = tci, where td is the cluster lifetime, typi- 
cally 100 Myr (Lada & Lada 2003). The final value of 0.1 pc'^ 
is the field star density as measured in the solar neighbourhood 
as well as in relatively old open clusters (Abt, 2009). The linear 
dependence adopted is based on the study of the dynamical evo- 
lution of the Pleiades over 1 Gyr by Converse & Stahler (201 1). 
They found that the central surface density (pC^) of the cluster 
falls linearly with time at a rate such that the initial density drops 
by tenfold in 200 Myr, and then falls also linearly but at a slower 
rate until 1 Gyr when the cluster is totally dispersed. They also 
show that the core radius is almost constant during the first 125 
Myr such that the linear dependence of the surface number den- 
sity (pc"^) transfers to the volume density (pc^^). We also note 
that in a similar way, Kaib & Quinn (2008) use this a priori time 
dependence for their sudy of the formation of the Oort cloud in 
an open cluster environment. 

The number of stellar encounters of miss distances d < dene 
can be tracked by the phase (cycle) 



1 {d < dene) = 



^ I ^enc " ^ 



(4). 



Stellar flybys of miss distances < dene, mean velocity Acrj y/n, 
and in the environment of varying star density n of eq.(3), occur 
successively at each integer value taken by this phase. 

We checked our algorithm in predicting the close stellar 
encounter rate F modeled by Proszkow & Adams (2009) in 
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their N-body simulations of stars in an embedded cluster Their 
power law F = Fq X (b/lOOOAUY corresponds to the number 
of close stellar encounters of miss distances < b per star per 
Myr. With our model, we estimated the time Tiqqoau required 
for the first stellar flyby at li < IQQQAU to occur by mak- 
ing cf> (d < lOOOAU) = J^'"""'" ijt^^^ dt = 1 with the con- 
stant stellar density n - 1000 stars pc^^ in eq (2). In a simi- 
lar way, we also estimated TsmAU for d < 500AU. With these 
two values Tufx^AU and Tsooau, we determined Fq = 0.4 flybys 
ib < 1000 AU) per star per Myr and y - 1.5 if cr = 1 km/s, and 
determined Fq = 1 flyby {b < 1000 AU) per star per Myr and 
y - 1.9 if cr = 5 km/s. These values of Fq and y are consistent 
with those in Proszkow & Adams (2009). As an additional test of 
our algorithm, we also predicted for the Sun that a close stellar 
encounter of 0.5 pc occurs once every 1.2 Myr in the present- 
day star field (n=constant=0.1 star pc~^ in eq (2)). This predic- 
tion is indeed the frequency determined by Garcia-Sanchez et al. 
(2001). 

3.2. Depletion of debris disks 

We turn now to our inital problem and apply Eqs. (2), (3), and 
(4) to estimate the fraction of planetesimals left in a debris disk 
after close stellar encounters have occurred for the first 100 Myrs 
of its lifetime. For this estimate, we need to determine the miss 
distance for which six stellar flybys occur during this period, 

l/fe„c dt - 6 with 
tci = 100 Myr We demonstrate below why six flybys are re- 
quired. We solve numerically this equation for dc„e and use its 
value in Figs|2]or[3]to estimate the fraction of stripped planetes- 
imals. This fraction is a lower limit, i.e. more planetesimals are 
ejected in reality, for two reasons. First, the miss distance of an 
encounter tracked by eq (4) is dene or smaller, and closer encoun- 
ters more severely deplete a disk. Second, stellar flybys at larger 
miss distances that are tracked hy (p (d < de„e) = 12, 18,24, ... 
in our algorithm are ignored in this work, although these flybys 
do slightly increase the fraction of stripped planetesimals. This 
ensures that our final conclusion is conservative but not overly. 
In practice, we solve numerically the equation <p (d < dene) - 6 
above for each stellar spectral type by using the corresponding 
fraction fp of Table|2]that set the star number density n in eq. (3) 
and hence in l/te„e. The fractions fp account for the fact that 
massive passing stars are rare, but destructive for a disk, and 
low-mass passing stars are common, but less harmful. 

We return to the choice (p (d < dene) = 6, which is based on 
the efficiency with which flybys strip a debris disk depending on 
the inclination i of the passing star orbit as well as the sense of 
circulation as discussed in Sect. 2. As already pointed out, a stel- 
lar flyby strips a debris disk according to the fractions given in 



Table 2. Stellar mass distribution of passing stars 



Table 3. Closest stellar encounters of the Sun and corresponding 
fractions of stripped planetesimals *^*' 



Spectral type m,, 


Mass range 

(Mo) 


Fraction /,, 


M8-M5 


0.10-0.21 


43.0 % 


M4-M0 


0.21 -0.47 


31.5 % 


K8-K0 


0.47 - 0.80 


12.4 % 


G8-G0-F0 


0.80- 1.70 


8.5% 


A8-A0 


1.70-3.20 


2.7% 


B8-B5 


3.20 - 6.50 


1.3% 



Spectral type nip 


nip/inr 


dene 

(AU) 


Fraction Sj 


M8-M5 


0.15 


180 





M4-M0 


0.25 


200 


2% 


K8-K0 


0.5 


230 


3% 


G8-G0-F0 


1.0 


500 





A8-A0 


2.0 


875 





B8-B5 


4. 


1250 






'*> duration=100 Myrs; m^ = IMq; no = 3000pc"^ cr = 5 km/s; 
<p{d < de„c) = 6. 



Figures|2]or[3]only if the passing star and disk motions are pro- 
grade and if / < 45° (two conditions). Assuming that the angu- 
lar momentum vector of the passing star is randomly distributed 
over the whole celestial sphere, the differential probability that 
this vector goes through a celestial band between inclinations / 
and ; + di is dp - jsin(i)di. By integrating between / - and 
7r/4, the probability corresponding to the two conditions above 
is p ~ 1/6. Hence, six flybys must occur to have one that is ef- 
ficient at stripping the planetesimals of a disk according to the 
fractions given in Figures |2]or[3] 

Finally, we combine the stripping effects over the six stellar 
spectral type ranges defined in Table |2] by computing the total 
fraction of the planetesimals left in the disk after 100 Myr, (1 - 
.Si) X (1 - ^2) X .... X (1 - .Sf,), where each 5, is the fraction of 
stripped planetesimals of Figures |2]or[3]for the appropriate ratio 
ntp/me. 

As an example of this calculation, we list in Table[3]the miss 
distance dene and stripping factor .?, for each spectral type of the 
stellar encounters undergone by the Sun for the first 100 Myr 
after its birth in a cluster characterized by cr = 5 km/s and 
the initial stellar density «o = 3000 pc^^ (Adams & Laughlin, 
2001). In these conditions, the total fraction of planetesimals left 
is < 95% after 100 Myrs in our model. The present day Kuiper 
Belt has only ~ 1 % of the mass that the minimun mass solar 
nebula hypothesis predicts, and, the hypothesis of the late heavy 
bombardement at 700 Myr has been put forth as an explaination 
(Morbidelli et al 2005). 

In Table H) we search the parameter space (nie, no) to find 
out the conditions in which a debris disk can be most severely 
depleted during the first 100 Myr of its lifetime. The frac- 
tions of planetesimals left after 100 Myr in a standard size disk 
(40 - 100 AU) and compact disk (10-40 AU) are estimated for 
the central star masses, me - 0.25, 0.5, 1.0 or 2.5 Mq, which 
cover the mass range of stars searched for debris disks in sur- 
veys. The cluster dispersion velocity adopted is cr = 5 km/s, 
as observed in 11 nearby open clusters and associations with 
ages between 5Myrs (Upper Sco) and 757 Myrs (Praesepe) by 
Madsen, Dravins & Lindegren (2002), and as measured in the 
N-body simulation of the dynamical evolution of embedded 
clusters from which open clusters emerge (Proszkow & Adams 
2009). In Table H) the initial star number density no is chosen 
to be 100, 1000, 3000, 10 000, 20 000, or 30 000 pc'^ for sev- 
eral reasons; most embedded clusters have stellar densities of 
~ 100 pc-^ (Carpenter 2000, PoiTas et al 2003, Lada & Lada 
2003), the Sun is thought to be born in a cluster of stellar den- 
sity ~ 3000 pc^^ (Adams & Laughlin 2001), and the densest 
and closest embedded cluster is the Orion nebula cluster with 
~ 20 000 pc-^ (Hillenbrand & Hartmann, 1998). We note also 
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Table 4. Fractions of planetesimals left in a debris disk after 
undergoing close stellar encounters for 100 Myr in expanding 
open clusters of various initial star densities. 



Central star 


Cluster**^ 


Fraction of planetesimals left 


trie (A^q) 


no (pc"') 










Disk 40-100 AU 


Disk 10 - 40 AU 






(%) 


(%) 


0.25 


100 


100 


100 


0.5 


" 


100 


100 


1.0 


" 


100 


100 


2.5 


" 


100 


100 


0.25 


1000 


100 


100 


0.5 


" 


100 


100 


1.0 


" 


100 


100 


2.5 


" 


100 


100 


0.25 


3000 


<65 


100 


0.5 


" 


<83 


100 


1.0 


" 


<95 


100 


2.5 


" 


100 


100 


0.25 


10 000 


<23 


100 


0.5 


" 


<49 


100 


1.0 


" 


<75 


100 


2.5 


" 


<94 


100 


0.25 
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that the central stellar density of the Arches cluster close to the 
Galactic center reaches ~ 10^ pc^^ (Portegies Zwart et al. 2007) 
but we have not included this extreme case in our study. 



4. Discussion 

According to our search of the parameter space {nic, no) in Table 
m severe depletion by close stellar encounters occurs only for 
a disk of standard size (40 - 100 AU) surrounding a star born 
in an embedded cluster with a high star-density uq greater than 
20 000 pc^^. In these conditions, fewer than 58% of the planetes- 
imals are left around an intermediate-mass star after 100 Myr, 
fewer than 37% around a solar-mass star, and fewer than only 
3% around a low-mass star In common low star-density em- 
bedded clusters where no is < 1000 pc^^, stellar flybys have 
a relatively small effect on disks. The turning point in Table |4] 
where disks start to lose their planetesimals in 100 Myr is the 
intermediate star density no ~ 3000 pc^^, which is thought to 
have prevailed in the birthplace of the Sun. In contrast, disks of 
compact size (10-40 AU) are almost insensitive to their stel- 
lar environment as seen in Table |4] These conclusions remain 
qualitatively the same regardless of whether the disk is dynam- 
ically excited at the start of the simulation or not. In agreement 
with these results, Spurzem et al. (2009) demonstrated that dis- 
ruptions of some wide-orbit planetary systems in an Orion-type 



■ 10 000 pc" 
20 000 pc" 
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0.1 0.2 0.5 1 2 

Mass of central star (Mg) 

Fig. 4. Fractions of planetesimals left in a debris disk after un- 
dergoing close stellar encounters for 100 Myr in open clusters of 
various initial star number densities. Fractions are the same as in 
Table|4]for the disk with inner and outer radii of 40 and 100 AU. 
The four bins correspond to the stellar spectral types M8-M5, 
M4-M0, KGF and A. 



cluster are expected on a timescale of a few 10^ yrs, leaving free- 
floating planets as relics. 

We can make a rough estimate of the fraction of stars born 
in low and high star-density environments by comparing the 
numbers of stars A^i and N2 that are, respectively, in the closest 
high star-density embedded cluster, the Orion Nebula Cluster at 
450 pc, and in all low star-density embedded clusters closer than 
this distance; these numbers are A^i - 2520 and N2 - 1 324 based 
on the catalog of Lada & Lada (2003). Thus, the corresponding 
fractions are N2/(Ni + N2) = 34% and A^i/(A^i + N2) = 66% 
for stars born in low and high star-density embedded clusters, 
respectively. If this rough estimate holds over the whole Milky 
Way, we can conclude that about two-thirds of the stars surveyed 
today originate from high star-density embedded clusters, thus 
their disks may have suffered severe stripping by close stellar en- 
counters according to Table |4] We note that embedded clusters 
with star-densities higher than 10 000 M^pc^^ survive the initial 
gas expulsion phase ("infant mortality phase"), and are stable 
against disruption by stellar evolution and encounters with gi- 
ant molecular clouds over an Hubble time (Lamers et al, 2005, 
Gieles 2009). 

As depicted in Figure [H or equivalently in Table H) the frac- 
tion of planetesimals left in a standard disk (40 - 100 AU) de- 
pends significantly on the mass of the central star nic- This is 
expected in our model in which flybys are more disruptive for 
disks around M-dwarfs encountering mostly more massive stars, 
while they are less damaging for disks around A-dwarfs encoun- 
tering mostly less massive stars. Consequently, in our model, 
disks surrounding M-dwarfs are expected to be more diflicult to 
detect since more depleted disks produce less observable dust. 
This trend in our model agrees with the observed decline in the 
number of debris disks found around lower mass stars in surveys 
(32% for A-dwarfs, 16% for FGK dwarfs, < 5% for M-dwarfs) 
as emphasized in Lestrade et al (2009). 

A limitation of our model is that the star number density 
is assumed to be uniform across the cluster Mass segregation 
is a well-known feature of star clusters, as can be seen in the 
Trapezium where the most massive stars are arranged in a cen- 
tral, compact core (Moeckel & Bonnell 2009). Mass segregation 
can be approximated in our model by a relative increase in the 
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high-mass star fraction fp given in Table |2] to simulate the con- 
dition at the center of the cluster For example, in the case of 
nic - 0.25 Mq and «(> = 20 000 pc"-', after increasing fp by a fac- 
tor of four for the most massive stars of this study {nip - 2.5M0) 
and decreasing fp accordingly for the lowest mass stars, the frac- 
tion of remaining planetesimals changes from < 3% in Table |4] 
to < 2%. For the same modifications of fp but for m^ = 2.5 M© 
and no - 20 000 pc"^, the fraction changes from < 58% in Table 
ID to < 45%. Consequently, this limitation of our model makes 
our estimate conservative and leaves unchanged our main con- 
clusion. 



5. Conclusion 

We have studied the depletion of planetesimals in a debris disk 
triggered by close stellar encounters in the environment of an 
expanding open cluster over its lifetime of 100 Myrs. We have 
found that depletion is significant only for an initial star-density 
of the embedded cluster of origin greater than 20 000 pc^^, as 
in the Orion Nebula Cluster, and in a disk of standard size (in- 
ner and outer radii : 40 and 100 AU). In these conditions, a de- 
bris disk loses >97% of its planetesimals around an M-dwarf, 
>63% around a solar-type star, and >42% around an A-dwarf 
in 100 Myr. This level of depletion could affect two-thirds of 
the stars searched in surveys since two-thirds of them are born 
in high star-density embedded clusters according to the catalog 
of Lada and Lada (2003). However, more compact disks (in- 
ner and outer radii : 10 and 40 AU) are much less sensitive to 
their environment. Unfortunately, debris disk sizes are not yet 
well enough known to decisively conclude whether the observed 
trend of fewer debris disks being detected around lower mass 
stars can be explained by the mechanism we have studied. 
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